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ABSTRACT 

Very Large Array and Parkes 64 m radiotelescope 21-cm observations of the 
starburst dwarf galaxy NGC 5253 reveal a multi-component non-axisymmetric 
H I distribution. The component associated with the stellar body shows evidence 
for a small amount of rotational support aligned with the major axis, in agree- 
ment with optically measured kinematics and consistent with the small galaxian 
mass. Approximately 20-30% of the H I emission is associated with a second 
component, an H I "plume" extending along the optical minor axis to the south- 
east. We consider outflow, inflow, and tidal origins for this feature. Outflow 
appears improbable, inflow is a possibility, and tidal debris is most consistent 
with the observations. Thus, kinematics of the H I that include this feature 
are not indicative of the dynamical mass or the local velocity dispersion of the 
cold gas. These H I observations also reveal a filamentary third component that 
includes an 800 pc diameter H I shell or bubble to the west of the nucleus, co- 
inciding with an Ha shell. The mass of H I in the shell may be as large as 
~4 X 10^ Mq. This large mass, coupled with the lack of expansion signatures 
in the neutral and ionized gas {v < 30 km s~^), suggests that this feature may 
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be an example of a starburst-blown bubble stalled by interaction with a massive 
neutral envelope. Many other H I kinematic features closely resemble those seen 
in Ha emission from the ionized gas, supporting the interpretation of neutral and 
ionized gas outflow at velocities of ~30 km s~^. Comparison between extinction 
estimates from the Balmer emission-line decrement and the H I column densities 
suggest a gas-to-dust ratio 2-3 times the Galactic value in this low-metallicity 
{Z = 1/4 Zq) galaxy. 

Subject headings: galaxies: individual (NGC 5253) — galaxies: ISM — galaxies: 
kinematics & dynamics — galaxies: starburst 



Introduction 



NGC 5253 is a remarkable starbursting dwarf galaxy with a prominent minor axis du st 
lane located in the nearby Centarus A/M 83 galaxy complex (IKarachentsev et al.l 120071 ) . 
At visual wavelengths, the impression of NGC 5253 is dominated by its network of ionized 
gas filarnents which extend a cross its minor axis and reach beyond the stel l ar distributio n 
hodgelEoGel : ICrahaml EoSlI : ICaldwell fc Philhpsl Il989l : iMarlowe et aP Il995l : iMartinI Eggsf ). 
X-ray observations reveal several small starburst-heated bubbles but no rn onolithic super- 
bubble or galactic wind (IStrickland fc StevensI Il999l : ISummers et al.ll2004l ). It harbors an 
extremely luminous, compact, obscured site of nuclear star for naation containing several star 
clusters with masse s from 10^ - 10^ M(7) haying ages < 10^ vrs (IGonzalez-Riestra et al.l 119871 : 



Calzetti et al.lll997l : IVanzi fc Sauvagell2004l : iMartm-Hernandez et al.ll2005[ )7 M uch of the nu- 



cleus is obscured by high levels of patchy dust extinction u p to Ay = 25 mag (lAitken et al. 



19821 : iMoorwood fc Glasslll982l : lAlonso-Herrero et al.l 120041 ). Radio- wave continuum and re- 
combination line studies have revealed a dominant luminous (10^ Lq), compact (r=l pc), 
young (< 2.3 Myr) star cluster that ap pears to be gravitationally bound and potentially 



desti n ed to evolve into a globular cluster (jBeck et al.lll996l : iTurner et al.lll998l : iGorjian et al. 



200ll : iMohan et al.ll200ll ). The centimeter- wave spectral energy distribution is consistent 



with an optically-thick thermal Bremsstrahlung sourc^ having densities exceeding 10^ 



cm 



and emission measures > 10*^ cm ^ pc powered by > 10^^ Lyman continuum photons s 



Caldwell &: Phillips! (119891 ) studied the star formation history of NGC 5253 and determined 
that the light was dominated by a young component with an age between 10^ and 10^ 



Turner. Beck, fc hJ |200d) term this object a "supernebula" while iKobulnickv fc JohnsonI |l999l ) 



use 



the name "ultra-dense H II regions" (UDHII) by analogy with the Galactic ultra-compact H II regions 
Wood fc Churchwelli (|l989h . 
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yrs. Recent studies have concluded that all of the UV-selected star clusters in NGC 5253 
have estimated ages <20 Myr, indicating that either the star formation episode began very 
recently or else that the timescale for d ynamical disruption of star clusters is very short 
( iTremonti et al.ll200ll : [Harris et al.ll2004j ). The overall gas- phase metallicity of NGC 5253 is 
12+log(0/if) ~ 8.16 (about 1/4 solar), but some of the central regions contain an excess of 



nitrogen, consistent with local "pollution" from W olf-Rayet star winds (IWalsh fc Roylll989 
Kobulnicky et al.l 119971 : iLopez-Sanchez et al.l 120071 ). 



When observed at low resolution, the distribution and kinematics of its neutral inter- 
stellar medium suggest that th e H I in NGC 5253 either rotat es about the minor axis or flows 
radi ally along the miri or axis (iKobulnicky &: Skillmarull995l ). iTurner. Beck. &: HurtI (119971 ) 
and lMeier et al.l (120021 ) report the ^^CO detection of molecular gas clouds having kinematics 
consistent with infall along the minor axis dust lane. These peculiar gas kinematics may 
hold clues to the mechanism that triggered the present burst of star formation. 

In this paper we present 21-cm aperture synthesis observations from the National Radio 
Astronomy Observatory's Very Large Arra^ (VLA) and the Parkes 64 m radio telescope to 
investigate the conditions in the interstellar medium that have produced such extreme star 
formation. NGC 5253 is an ideal environment to seek the causes of extreme star formation 
because the burst is so young and, presumably, feedback from the energetic star clusters has 
not yet had a dramatic impact on the surrounding ISM. 

The distance to NGC 5253 is s omewhat controve rsial and has been ra easured in many 
studies including ISaha et all Jl995h (4 1 ±0.5 Mpc). Icibson et all J200oh (3.3± 0.3 Mpc), 
Freedman et~aD J200lh (3^25 ±0.2 Mpcl.lKarachentsev et al.N2002h (3.9±0.5 Mpcl. lThim et al 
(I2OO3I) (4.0±a3 Mpc), and ISakai et all (|2004h (3.8 ±0.2 Mpc). The distances based on 



Ceph eid variables strongly depend on the sample select ion (IGibson et al.ll2000l : iThim et al. 



2OO3I ) and so we prefer the TRGB distance of 3.8 Mpc (jSakai et al.ll2004l ) which implies an 
angular scale of 18.4 pc arcsec"^. 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc. (AURA) under 
cooperative agreement with the National Science Foundation. 
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2. Observations 
2.1. VLA 

NGC 5253 was observed in tlie 21-cm line of neutral hydrogen using the VLA for 
approximately 6 hours in the BnA configuration on 1995 September 25, 3 hours in the CnB 
configuration on 1996 January 28, and 2 hours in the DnC configuration on 1996 June 
10. The 2IF correlator configuration with Hanning smoothing yielded a total bandwidth 
of 1.56 MHz (327 km s~^) and 128 12.2 kHz (2.58 km s~^) spectral channels. Intermittent 
observations of the compact radio source 1316-336 were used to phase calibrate the data, and 
observations of 3C286 were used for flux calibration. The data were reduced with NRAO's 
AIPS package following standard procedures for editing bad time ranges and performing 
phase, flux, and bandpass calibrations. The radio continuum flux was subtracted in the 
UV plane using offline channels to produce a continuum-free 21-cm line UV dataset. Data 
from each configuration were calibrated separately, and then the UV data were combined 
after subtracting continuum. These combined data were then used to map the emission 
from the 1420 MHz transition of neutral hydrogen with greatly enhanced angular resolution 
and sensitivity compared to the 45-minute DnC configuration observations presented by 



Kobulnickv fc SkillmanI fll995[ ). 



The Uy data were mapp ed to the image plane at two different resolutions and CLEANed 
( IClarkI Il980l : iHogbonj Il974j ) with resulting clean synthesized beamsizes of 9.0" x 7.6" and 
17"xl5" to emphasize both the small scale structures and diffuse H I features, respectively. 
The RMS noise levels in these maps are 1.2 mJy beam~^ and 1.7 mJy beam~^, respectively, 
slightly larger tha n the theoretical noise o f 1-0 m Jy. Because the angular extent of H I in 
NGC 5253 (~ 2'; iKobulnickv fc SkillmanI ffl995h ) is much smaller than the VLA primary 



beam size at 21 cm, no primary beam correction was applied. The low-resolution cube was 
used as a template in which regions having emission exceeding 2a in at least two adjacent 
channels were marked. The low resolution cube was then used to blank the higher resolution 
cube on a channel-by-channel basis outside the areas of probable emission. This approach 
preserves low-level 1-2 a emission peaks, many of them real, in the final high-resolution data 
cube used for subsequent analysis. 



2.2. Parkes 

NGC 5253 was observed with the Parkes 64 m radio telescope in the 21-cm line of 
neutral hydrogen on 1997 February 12 using a total bandwidth of 8 MHz (1600 km s~^) and 
2048 channels of 0.0039 MHz (0.82 km s~^) centered at Vq = 250 km s~^. The primary beam 
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angular diameter is 16'. We interleaved on-band spectra with frequency-switched off-band 
spectra centered at Vq = 880 km s~^ to obtain emission-free H I profiles. A 9-pointing 
(3x3) grid with 14' spacings centered on NGC 5253 was used to search for extended emission 
surrounding NGC 5253, but only the central position showed a strong 21-cm detection. The 
data for left and right circular polarizations were reduced by subtracting the frequency- 
switched reference spectrum from the source and fitting a low-order polynomial to remove 
any residual baseline. The system temperature averaged 33 K and the mean rms noise in 
the resulting spectra is 0.06 K per channel. We observed the Galactic H II region S8 as a 
flu x calibration st a ndard and scaled the observed temperatures to match the 76 K reported 



bv lKalberla et all fll982f ) 



3. Analysis of the New H I Observations 

3.1. Total H I and Line Profiles 

Figure [1] shows the 21-cm emission line profile (top panel) of NGC 5253 from the Parkes 
64 m telescope (dotted line), the VLA D configuration short-baseline data (solid line) and 
the VLA B+C configuration high-resolution data cube (dashed line) as described above. 
The lower panel shows the 21-cm absorption line profile from the longest baselines of the 
VLA B configuration against the co mpact thermal rad io continuum source associated with 



the super star cluster in the nucleus (ITurner et al.lll998l ). The integrated single-dish H I flux 



is 43.9 Jy km s ^, equivalent to a neutral hydrogen mass of 1.42 x 10^ Mq at a distance of 
3.8 Mpc, assuming. 



HI 



2.36 X 10^[D(Mpc)]2 J S{Jy)5v{km s'^). (1) 



The VLA D configuration data recovers 77% of this total, implying the presence of some 
3 X 10^ M© of neutral hydrogen distributed smoothly on scales larger than ~ 7', the largest 



struc tures visible in 2 hours to the VLA at 21 cm in the D configuration ( jUlvestad et al. 



20061 ) . The high-resolution data cube contains 42% of the single-dish flux, or about 6.0 x 10^ 
Mq. Therefore, H I column densities measured from this cube will be lower limits to the 
true values which may be larger by a factor of order < 2. 

The 21-cm absorption spectrum in the lower panel of Figure [1] shows a broad absorption 
line spanning at least 80 km s~^ with a maximum optical depth of r = 0.55 near Vq = 409 
km s~^. The velocity of maximum optical depth is redshifted ~ 10 — 20 km s~^ compared 
to the 21-cm emission profile peak at Vq ~ 395 km s~^. This is consistent with the central 
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radio supernebula lying behind the bulk of the cool atomic gas and the latter having a net 
radial vel ocity toward the n ucleus, as su ggested from the kinemat ics of the cold molecular 
clouds bv iMeier et all J2002h . Adopting Jpickev fc Lockman|[l99oh . 



AT, 



H 



T Av Ts 
5.2 X 10-19' 



(2) 



where r is the 21-cm optical depth, Tg is the spin temperature and Av is the velocity full 
width at half maximum of the 21-cm absorption line, we can estimate the column density, 
Nh, of neutral hydrogen in front of the compact nuclear radio supernebula. Because the solid 
angle probed by the VLA absorption and emission spectra are vastly different, we do not 
have an estimate of the sp in temperature for this gas, so we adopt a lower limit of > 50 K 
(jPickey fc Lockmanlll990l ). The overall width of the feature is consistent with an abundance 
of cold H I clouds in front of the nucleus over a range of velocities Vq = 370 — 435 km s~^. 
Therefore, we adopt a lower limit of Av > 10 km s~^. The resulting H I column densities 
are Nh > 5 x 10^° cm"^ at this location. 



3.2. The H I Distribution at High Resolution 



Figure [2] shows the integrated H I column density (contours) overlaid on a 6450 A optical 
continuum image of NGC 5253 (greyscale). The inset at lower left depicts the 9.0" x 7.6" 
synthesized beam. Contours show 21-cm line fluxes of 0.04, 0.08, 0.12, 0.16, 0.24, 0.36, and 
0.40 Jy km s-^ corresponding to beam-averaged H I column densities of 6, 12, 18, 24, 36, 
52, and 64 xlO^° cm"^. A cross marks the location of the dominant embedded star cluster 



and thermal radio supernebula (ITurner et al.lll998l ) where H I appears in absorption. The 



H I column density peaks near 6.4x10 cm at ~15" (270 pc) south-southeast of radio 
superne bula and near th e dust lane which coincides with the position of several molecular 
clouds jMeier et al.ll2002h . 



As noted by iKobulnicky &: SkillmanI (119951 ) . the H I distribution in NGC 5253 does not 
closely follow the optical morphology. Rather, the H I is extended along the minor axis to 
the southeast and northwest. Note that at this higher resolution, it is clear that the extended 
gas is not symmetrically placed along the minor axis. The extension to the southeast is a 
much more prominent feature, and we designate this as the H I "plume". These new data 
also reveal arc-like extensions of H I to the west of the nucleus, suggestive of a partial shell 
or bubble of neutral hydrogen. There are also extensions of H I to the NE and SW that 
generally follow the optical major axis. 
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Figure [3] depicts the integrated H I column density (contours) overlaid on an Ha image 
(greyscale). Here we note that the morphology of the H I shell to the west corresponds 
closely to the morphology of the Ha shell at position angle of about —70° as seen from the 
central radio supernebula. The morphology of the H I to the east appears to mimic the Ha 
emission in the inner parts of the galaxy where the H I peak column densities coincide with 
the Ha peaks. Also in the east, similar to in the west, there is an Ha extension that appears 
to be outlined in H I (at a position angle of about 100°). However, further from the center, 
this correlation breaks down. Specifically, the bright H I plume, at a position angle of about 
140° does not have a corresponding strong Ha emission feature. 

Figure H] depicts the integrated H I column density (contours and greyscale) at a lower 
resolution of 16" useful for identifying lower column density features. A few additional 
features appear in this figure. The H I extension to the west now appears more clearly and 
shows the same open shell geometry as the Ha. At this lower resolution, it is now possible 
to see that there is H I at a large radius from the center but on the same position angle 
as the center of the wes tern open shell. This loo ks reminiscent of simulations of bubbles 



that achieve break-out (iMac Low fc Ferraral Il999l ). Also, at lower resolution, the eastern 



extension looks more like the same open shell morphology and also follows the Ha. 

In summary, the integrated H I maps reveal three neutral gas components that will 
become more clear in the velocity channel maps below: 1) extensions to the NE and SW 
that follow the optical major axis, 2) an H I plume to the SE that extends well beyond the 
optical extent of the galaxy, and 3) open shells in the west and east, and 



3.3. The H I Kinematics 

Figure [5] shows 21-cm maps of NGC 5253 in every 2.58 km s~^ velocity channel contain- 
ing emission from Vq = 459 to Vq = 358 km s~^. The greyscale range shows beam-averaged 
column densities of 1.0x10^° cm~^ (white), equivalent to 2x the rms noise, to 5.0x10^° cm~^ 
(black). The cross marks the location of the central super star cluster and radio supernebula. 
Ellipses indicate the positions of the mo lecular clouds A th rough E (from left to right, in 



decreasing Right Ascension) identified by iMeier et al.l (120021 ) . In Figure [5l the southeast H I 



plume is the dominant feature at higher velocities (450 - 420 km s ^), while the H I aligned 
with the stellar distribution is found primarily at velocities between 420 and 380 km s~^. 

The kinematics of individual H I features relative to the ionized gas can be better seen 
in Figure [6] which shows H I column density channel maps (contours) overlaid on an Ha 
image (greyscale). Every second 2.58 km s^^ channel is plotted. The contours show fluxes 
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of 2.4, 3.6, 6.0, and 8.4 mJy beam~^, equivalent to beam-averaged H I column densities of 
1.0, 1.5, 2.5, and 3.5 xlO^° cm~^. Again, it can be seen that the highest velocity gas at 
Vq = 420 — 455 km s~^ (upper 2 rows of panels) is located to the southeast of the nucleus. 
Gas at intermediate velocities Vq = 400 — 420 km s~^ near the systemic velocity of Vq ~ 400 
km s~^ has a circularly symmetric morphology similar to that of the Ha. At these velocities 
there are also lower column density H I extensions along the major axis of the galaxy to the 
northeast, southwest and also in the direction of the Ha shell to the west. At the lowest 
velocities, Vq < 400 km s~^, the H I morphology becomes elongated along the galaxy's major 
axis. 

Figures [3] and [6] also show H I extensions to the northeast and southwest along the 
galaxy's major axis. The H I protrusion to the northeast appears at velocities Vq = 415 — 380 
km s~^ and does not unambiguously form a coherent shell or bubble. The H I extension 
to the southwest appears over a similar velocity range with more gas closer to the systemic 
velocity. The spatial and kinematic symmetry of these extensions along the major axis is 
consistent with low amplitude rotation, in the canonical manner, about the galaxy's minor 
axis, but the pattern is, by no means, obvious. 

Figure [7] shows a 21-cm intensity-weighted velocity map ( 1st moment; contours) overlaid 
on a stellar continuum 6450 A image (greyscale). Isovelocity contours are labeled in km s^^. 
This figure shows a steep kinematic velocity gradient along the southeast H I plumejf] The 
contours show the largest radial velocities of Vq ~ 450 km s~^ at a distance V to the south- 
southeast of the nucleus. Radial velocities approach the systemic velocity of Vq ~ 395 km 
s~^ near the nucleus. Along the major axis to the northeast and southwest of the nucleus, 
the velocity field is disordered and shows a small velocity spread, ranging between Vq ~ 385 
in the northeast and Vq ~ 405 km s~^ in the southwest. 



Kobulnicky fc SkillmanI (119951 ) concluded, on the basis of low-resolution (45") H I maps, 
that the neutral atomic medium in NGC 5253 rotates about the major axis, as in polar ring 
galaxies. However, these new high-resolution H I data reveal a strong asymmetry in the H I 
distribution and velocity field, making the rotation scenario improbable. These new data 
are more consistent with an either an inflow or outflow of neutral gas along the minor axis 



from t he south and east o f the nucleus. Inflow, was co nsid ered bv iKobulnickv fc Skillman 



(I1995I ) and propounded by lTurner. Beck, fc HurtI (119971 ) and lMeier et al.l (120021 ) on the basis 



of ^"^CO observations. 



•^The gradient is ne arly parallel to the projected direction toward the much larger group member M 83, and 
van den Berghl |l98d) suggested a recent interaction between the two. However, the line-of-sight separation 



exceeds 1 Mpc given the current dis tance determinations of 3.8±0.2 Mpc for NGC 5253 and 5.16±0.41 Mpc 
for M 83 (jKarachentsev et alJl2007l l. effectively precluding any recent interaction. 
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Figure [8] shows the 21-cm intensity-weighted velocity (1st moment) map as in Figure [7] 
(contours) overlaid on an H I total intensity map (0th moment; greyscale). The scale bar 
shows H I surface brightness in units of Jy beam~^ m s~^. A peak flux of 400 Jy beam~^ m s~^ 
corresponds to H I column densities of 6.4x10^^ cm~^. The peak H I column densities occur 
where the velocity gradient is the steepest. The majority of the H I by mass is located to 
the south-southeast and in the circumnuclear region. H I extensions to the west, southwest, 
and northeast have comparatively low total column densities and low line-of-sight velocity 
dispersions, consistent with gas associated with the main body of the galaxy. 

Position-velocity diagrams allow another perspective on the observations. Figure [9] 
shows a position-velocity diagram of the H I emission taken along the optical major axis 
at position angle 50°. Although the diagram is quite complex, two main features stand 
out. First, there is a large velocity range in the central parts. Here one can appreciate the 
asymmetry in the field. Redshifted H I gas extends out to 50 km s~^ from the systemic 
velocity, while blueshifted gas extends out only 30 km s~^. The second notable component 
is the gas at roughly the systemic velocity showing a slow drift in central velocity across the 
spatial distribution with an amplitude of roughly 20 km s~^, consistent with that expected 
for such a low-mass galaxy. 

Figure [To] shows a position- velocity diagram of the H I emission taken along the optical 
minor axis at position angle 140°. Here we again see a large range in velocities over a relatively 
small spatial extent. It is interesting that in this projection, the gas appears continuous in 
velocity from the extreme of the redshifted gas to the extreme of the blueshifted gas. 



4. Interpretation of Major H I Components in NGC 5253 

Here we will try to draw a coherent picture of the three major components of the H I 
distribution: the H I associated with the stellar distribution, the southeast H I plume, and 
the open shells to the west and east. 



4.1. HI Along the Major Axis 



It is interesting to compare these new H I observations of NGC 5253 to those of NGC 625 
( ICannon et al.ll2004l ). another dwarf starburst galaxy at a similar distance. ICote et al.l (120001 ) 
showed that NGC 6 25 also has the dominant H I velocity gradient along the minor axis. 
Cannon et al.l (120041 ) showed that this large velocity gradient was due to outflowing H I gas 
from the central starburst. This outflowing gas is also detected in the warm and hot phases 
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of the ISM flCannon et al.l 120051 ) . 



This similarity with NGC 625 motivates a second look at the H I component in NGC 5253 
that is aligned with the optical galaxy. In NGC 625, the signature of rotation along the ma- 
jor axis was swamped by the presence of the strong H I outflow. However, the position 
velocity diagrams revealed the low amplitude rotation curve present in the H I. It is proba- 
ble that we are seeing the same effect in NGC 5253. Although the H I in NGC 5253 is quite 
disturbed, the overall trend for higher velocities in the southwest and lower velocities in the 
northeast suggest a patter n of rotation. This is also seen in the stellar velocities meas ured by 
Caldwell fc PhiUipsI Jl989h (their Figure 10) and the Ha major-axis velocity map of iMartinI 



(jl998l ) (her Figure 9). It is not straightforward to calculate the expected amplitude of the 
rotation curve for NGC 5253; comparing to the luminosity requires a significant correction 
for the small M/L ratio expected for a starburst galaxy. Nonetheless, after applying a small 
(< 20%) correction for inclination, the amplitude of ~ 10 km s~^ is not out of line with what 
could be expected fo r NGC 5253. Note that the H I study of the similar dwarf starburst 
galaxy NGC 1569 by lStil fc Israeli (120021 ) also found only small evidence of normal rotation 
in the inner parts of NGC 1569, but that the H I along the major axis in the outer parts 
confirmed this rotation signature. 



4.2. The H I Plume 



These new high resolution H I observations preclude rotation along the minor axis as a 
viable interpretation for the southeast plume. Either inflow, outflow, o r some kind of tidal in 



teract ion remnant feature seem more likely. Based on CO observations, [Turner. Beck, fc Hurt 



( 119971 ) favored infall of metal poor gas. The radial velocities of the CO clouds ~200 pc east 
of the nucleus are slightly higher than the systemic velocity, consistent with inflow or pos- 
sibly rotation about the major axis within the innder few hundred pc. The H I velocities 
at this location could also be consistent with major-axis rotation within the nuclear regions. 
Many galaxies are known to have counter-rotating or decoupled dynamical components in 
the inner regions, and close grvitational interactions are often invoked as the causes of such 
kinematics. 

That an accretion event may have triggered the current starburst ha s been suggested sev- 
eral ti mes Jvan den Berghlll98ol : lGrahanJll98ll : [Caldwell fc PhiUipJl989] l iTurner. Beck, fc Hurt 
( 119971 ) identify three points in favor of the infall hypothesis. First, the weak CO emission 
associated with a strong starburst favors low metallicity gas. Second, the non-central posi- 
tion of the CO suggests a transient feature. Third, t he redshift of the C O relative to the 
galaxian systemic velocity could be indicative of inflow. iMeier et al.l (120021 ) further note that 
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since the dust lane that is associated with the CO and H I emission in the southeast plume 
appears to be in front of the galaxy, the geometry (redshift) favors infall. 

A close look at the relative locations of the CO clouds, dust lane, and H I peaks along 
the southeast minor axis reveals, in fact, two different gaseous components in this direction. 
These two components are resolved both spatially an d kinematically by th e new H I observa- 
tions, and a comparison with the extinction maps of ICalzetti et al.l (119971 ) shows that these 
are intrinsically different components — not two maxima of the same feature. 

Figure [11] shows a gre yscale representation of the dust reddening (Ay) computed from 
the Ha/ Hf3 ratio maps of ICalzetti et al.l (119971 ) assuming Ay = 3.1E{B — V). The inferred 
extinctions range from Ay{Ha/ Hf3) near zero to Av{Ha/Hj3) ~ 2. The cross and ellipses 
show the locations of the central thermal radio supernebula and the CO clouds, as in previous 
figures. Contours in Figure [TT] show the extinction derived from the ne utral hydrogen column 
densities assuming a Galactic ratio of gas to dust (IBohlin et al.l 1 19781 ) . 



{Nhi + 2Nh,] 
E{B - V) 



5.8 X lO^^atoms cm ^ mag ^ 



(3) 



For a standard Galactic ratio of total to selection extinction, R = Ay / E{B — V)= 3.1, this 
is equivalent to. 



3.1 



{Nhi + 2Nh,) 



5.8 X lO^^atoms 



cm 



mag' 



A 



V 



(4) 



Here, we have ignored, for the moment, any contribution from molecular gas, and we 
remind the reader that the H I column densities inferred from the high-resolution interferom- 
eter maps may be underestimates of the true column densities by up to a factor of two, owing 
to the insensitivity of the VLA to emission on large angular scales. Hence, the estimated 
H I column densities and extinctions are lower limits. The contours in Figure [TT] correspond 
to Ay{HI) =4, 3.5, 3, 2.5, 2, 1.5, 1 and 0.5 mag. 

Figure [11] shows that the extended region of high extinction indicated by the Balmer 
decrement coincides with a local maximum in the neutral hydrogen column density and the 
locus of the positions of the CO detections. The ratio Ay{HI) / Ay[Ha / H [3) in this region 
~ 3.5/2 = 1.75, consistent with the expectation that some of the H I lies behind the region 
of ionized gas, and with the expectation that some of the diffuse H I may be unseen in the 
interferometer maps. More precise estimates are no t possible because of the large disparity 
in beamsize between the VLA H I data (8") and the ICalzetti et al.l (119971 ) HST maps (0.3"). 



Nevertheless, the correlation between both extinction measures at this location is consistent 
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with the majority of the H I (and molecular gas?) lying in the foreground of the ionized 
medium. The high H I column density in this feature probably continues into the region 
surrounding the central radio supernebula and super star cluster, but measurements are 
unreliable at this location owing to the strong H I absoprtion. 

By contrast, the maximum H I column density, corresponding to Ay {HI) = 4 mag, is 
located in the H I plume ~15" to the southwest of the region with the largest Av{Ha/ H (3). 
At this location there is no pronounced peak in Av{Ha/ HjS) and there are no CO detections. 
Figure shows that the gas responsible for this primary peak in N(H I) occurs at velocities 
Vq > 410 km s~^ and is kinematically and spatially distinct from the neutral gas at the 
location of the molecular clouds. Thus, what we are calling the H I plume appears to have no 
physical relation to the region of high extinction in the dust lane, with the possible exception 
of a tenuous connection at velocities Vq < 225 km s~^. 

The very high ratio of H I column density to extinction in the plume implies that ei- 
ther the plume gas is located on the far side of NGC 5253, or that the gas is extremely 
metal-poor and has a very large gas-to-dust ratio. The metal-poor ISM in dwarf galax- 



ies d oes exhibit gas-to-dust ratios 2-5 times higher than in the Milky Way (jPraine et al. 
20071 ). so equation H] could require a correction factor of this order. However, the mean 
ratio Av{HI)/Av{Ha/H(3) = 3.75/0.5 = 7.5 in the vicinity of the plume compared to 
Av{HI)/Av{Ha/H(3) = 3.25 near the dust lane suggests that the plume and the dust lane 
are distinct entities. Even if we allow that, statistically, half of the neutral medium may lie 
behind the star-forming ionized gas, the inferred extinctions are factors of 2-3 larger from 
the H I column density than from the Balmer decrement. The real ratio may be even larger 
because, as described above, the H I column density estimates should be regarded as lower 
limits and a molecular component is not included. 

The true nature of this low-extinction, high-column-density, redshifted H I plume is 
difficult to pin down. One possible explanation is that the H I pl ume may be outflo wing 
gas such a s commonly observe d in dwarf starbursts like NGC 1705 (IMeurer et al.lll998l ) and 
NGC 625 JCannon et al.ll2004[ ). The main difference between NGC 5253 and NGC 625 is the 
asymmetry of the H I distribution seen here. The symmetrical H I distribution in NGC 625 
was taken as strong evidence against infall of gas, but the H I plume in NGC 5253 does 
not show a similar symmetry. Figure [10] shows that there is continuity from the strongly 
redshifted gas (~ -|-50 km s~^) through to blueshifted gas (~ —30 km s~^), which would be 
the signature of s omewhat asymmetri c outflow. The "spur" of outflowing gas in NGC 1705 
is one-sided, and iMeurer et al.l (119981 ) note that this is to be expecte d when the source of 
the outflow is displaced relative to the midplane (IMac Low et al.lll989l ). However, the plume 
in NGC 5253 contains 20-30% of the H I seen in the interferometer maps, depending on 
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the spatial and velocity boundaries adopted. Such a fractionally massive outflow of neutral 
gas over a fairly small solid angle without any ionized gas counterpart makes the outflow 
interpretation less plausible in this galaxy. 



Radial inflow, suggested by iMeier et al.l (120021 ) on the basis of the kinematics of CO 



clouds within the minor axis dust lane, is another possibility. The geometry of the isovelocity 
contours in the redshifted H I plume can help distinguish whether the inflow could be on 
the near or far side of NGC 5253. The most redshifted contours (455 km s~^) lie farthest 
from the galaxy and are ~50 km s~^ larger than the systemic velocity of 400 km s~^. An 
infalling stream is expected to accelerate as it approaches the galaxy. Thus, if the plume 
were infalling gas on the near side, the isovelocity contours closer to the galaxy ought to 
exceed 455 km s~^. Instead, we observe that they approach the systemic velocity. Inflow 
from the far side of NGC 5253 might, however, produce the observed signature. Such a 
background geometry would simultaneously explain the lack of extinction seen toward the 
high column density peaks in the plume. 

Finally, and perhaps most plausibly, the H I plume may be a kinematically distinct 
remnant from a recent gravitational interaction with another member of the Centarus A 
subgroup of the Centarus A/M 83 complex. It is difficult to prove or rule out this possibility. 
In this scenario, the H I asymmetry is easily understood as a chance projection against the 
main body of NGC 5253, and there is no need for the tidal tail to be symmetrically placed 
in either space or velocity. The galaxy density within the Centarus A subgroup is high, and 
interactions are to be expected. Many of these systems are dynamically disturbed systems, 
suggestive of a history fraught with close encounters. A perturbing, but not destructive, 
close encounter would simultaneously provide an explanation for the H I plume as a tidal 
feature and provide a cause for the recent burst of star formation seen throughout the galaxy, 
most exceptionally in the nucleus. 

Karachentsev et al. (2007) have provided a three dimensional picture of the Cen A/M83 
complex, and NGC 5253 is located in the periphery. Using radial velocities and accurately 
measured distances, Karachentsev et al. (2007) have identified NGC 5253 with the Cen A 
sub-group based on a positive tidal index. This is because, even though NGC 5253 is within 
2 degrees of M 83 in projection on the sky, it is 1.4 Mpc closer than M83 and nearly identical 
in distance to Cen A. Note, however, that NGC 5253 is separately by roughly 12 degrees 
in projection from Cen A, corresponding to roughly 800 kpc. Thus NGC 5253 is currently 
far from the strong gravitational influence of either of the largest galaxies in the complex. 
Despite being a member of a very populous complex, NGC 5253 appears to be relatively 
isolated. In fact, the next nearest known galaxy, ESO 383-087, appears to be over 400 
kpc away. It is not clear which member of the group might have interacted recently with 
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NGC 5253. 



In some ways, the H I plume in NGC 5253 is reminiscent of some of the ex tended 
filaments of gas surrounding IC 10 ( Shostak fc Skillman 19891 ;^ Wilcots fc Miller 1998 ) or the 
giant H I cloud in Virgo, H I 1225+0146 jMcMahon et aljll99oh . In IC 10, H I "plumes" 
and "streamers" have velocity gradients that show little relation to the overall rotation of 
the galaxy. The true nature of the extended H I features in IC 10 is also not yet understood, 
but interac tion wi t h the se H I features is generally thought to have provoked the present 
starburst. iTaylorl (119971 ) finds that starbursting dwarf galaxies like NGC 5253 are found 
preferentially to have close companions, so the interaction hypothesis is a reasonable one. 

One very important point is that the large range of velocities observed in the H I in 
NGC 5253 is associated with the southeast plume. It is common practice in the literature 
to apply the FWHM of t he gas velocity dis tribution as indicative of internal dynamics t o 
estimate dynamical mass (IHarris et al.ll2004l ) or a stellar crossing time (jCalzetti et al.lll997l ). 
Regardless of whether the large-scale velocity gradient seen in the H I plume is due to infall, 
outflow, or tidal debris, the resulting large velocity dispersion should not be used as an 
internal dynamical indicator. 



4.3. The H I Shells and Bubbles 



The main question concerning the open shells on the west and the east is whether those 
shells represent events centered at the currently indicated positions or if they are also the 
result of centrally energized outflows that give rise to shell geometries as they emerge from the 
galaxy. The centers of the western and eastern open shells lie well outside of the main stellar 
distribution and far from the current, centrally concen trated star formation. Comparison 
with simula tions of outflows (IMac Low &: Ferraralll999l ) and real outflows observed in Ha 
and X-rays (jStrickland et al.l 120041 ) show that open-shell morphologies can be created where 
the centers of the shells are far from the source of the energy. Is it possible that this is what 
we are seeing in NGC 5253? 

How likely is it that star formation events capable of creating large H I shells could 
have occurred on the periphery of the present stellar distribution? Ou r view of the spatially 



resolved star formation history of NGC 5253 has evolved with time (Ivan den Berghl Il980 



Caldwell fc Phillipslll989l: 



most recent work by 



Calzetti et al. 



Harris et al 



1997: Tremonti et al. 2001; Harris et al. 



20041 ). The 



(120041 ) emphasizes the very youn g ages (< 20 Myr) fo r 
all of the clusters that they s tudied (in agreement with the work by lTremonti et al.ll200ll ). 
However, [Harris et al. required detection in the UV for inclusion, so their sample 
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should be strongly biased toward the youngest clusters (only the brightest of the older 
clusters would be detected - see their Figure 15). The earlier study by lCalzetti et al.l (119971 ) 
derived ages of order 50 Myr for some of the clusters, the c o lors o f the core required starburst 
activity within the last 100 Myr, and ICaldwell fc Phillips! (119891 ) deduced a range in cluster 
ages with an upper limit of 10^ years. From the above, it is not clear whether strong star 
formation events in the periphery of its stellar distribution in the recent past are consistent 
with the star formation history of NGC 5253. A star formation history derived from the 
field stars, while challenging owing to the distance and the presence of differential extinction, 
would be a valuable study. 

Figure [3] depicts the integrated H I column density (contours) overlaid on an Ha im- 
age (greyscale). The morphology of the H I shell to the west corresponds closely to the 
morphology of the Ha shell at position angle of about —70° as seen from the central radio 
supernebula. The shell has an angu l ar dia meter of ~45" or about 830 pc. This Ha feature 
lies along echelle slit #2 of iMartiru (119981 ) who noted the abser ice of k i nemat ic splitting in 
Ha that would normally demarcate an expanding ionized shell. iMartinI (119981 ) characterizes 
this as a quiescent shell. The mass of neutral hydrogen in this H I feature is ~ 4.8 x 10^ M© 
as measured by summing emission in the high-resolution data cube from this feature. Such a 
significant mass of neutral gas could be re sponsible for stalling the shell's expansion and pre- 



venting "blowout" . ICalzetti et al.l (120041 ) show emission-line ratio maps for the ionized gas 



in this vicinity, and they concluded on the basis of high [S II] /Ha ratios that the numerous 
shells and radial filaments to the south and west are shock ionized, consistent with starburst 
powered shells. However, a softer radiation field, characteristic of the diffuse i onized gas at 
large distances from star forming regions, could also explain these line ratios (IMartiru 119971 : 



Hunter k Gallaghei]ll990h . 



In Figure it can be seen that the H I coincident with the western Ha shell spans a 
range of velocities Vq = 380 — 410 km s~^. However, the fragmentary morphology of the 
neutral gas at these velocities does not provide clear evidence of a single coherent structure. 
The velocity range (30 km s~^) of the H I near this Ha shell is comparable to that found in 
the ionized gas, consistent ^vith sl ow or stalled expansion. The lack of X-ray emission inside 
this bubble ([Summers et al.ll2004l ) is also consistent with it being an older, cooler, quiescent 
remnant from bygone era of star formation. If we take the mass of the H I bubble to be 
4.8 X 10^ M© and the radial expansion speed to be 20 km s~^ the kinetic energy contained 
in this structure is 1.7 x 10^^ erg. The implied expansion timescale for this ~1 kpc diameter 
structure is <40 Myr, making it unlikely that the current 3 Myr old burst produced this 
bubb le. It is probably a pro duct of the older episodes of star formation from the previous 50 
Myr. iMarlowe et al.l (119951 ) estimate the mechanical energy injection rate for NGC 5253 to 
be 3.6 X 10^° erg s~^, or about 5 x 10^^ erg over the last 5 Myr. Even though star formation 
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and its effects are distributed throughout NGC 5253, this energy injection rate is more than 
sufficient to power the observed H I shell, presuming tha t such an inject i on ra te has persisted 
for at least 50 Myr — the ages of the clusters studied by lCalzetti et al.l (Il997l ). 



5. Comparison of the H I with Other Gas Phases 
5.1. Comparison with the Molecular Gas 



Figure [5] shows the positions of CO detections identified by lMeier et al.l (120021 ) as molec- 
ular clouds A through E (from left to right, in decreasing Right Ascension). Molecular cloud 
D is nearly coincident with the central radio supernebula (marked by a cross). H I appears 
in absorption at this location, so it is not possible to identify molecular cloud D with a 
particular 21-cm feature. Molecular clouds B {vlsr = 421 km s~^) and C {vlsr = 419 km 
s"^) have probable H I counterparts seen in Figure [5] between Vq = 430 km s~^ and f© = 405 
km s~^0 Molecular cloud A {vlsr = 363 km s~^) has no apparent 21-cm counterpart in a 
similar velocity range. However, H I emission is present at that location over the velocity 
range Vq = 418 — 395 km s~^. Molecular cloud E {vlsr. = 424 km s~^) has no corresponding 
H I counterpart in a sir ailar velocity range . Detectable 21-cm emission is present at velocities 
Vq = 410 — 380 km s~^. iMeier et al.l (120021 ) note that clouds D and E are tentative detections 
at the 3-5 a level of significance, and they may not be real. The lack of an H I counterpart 
at the velocity of cloud E is consistent with it being noise in the CO data. 



The peak H I column density of 6.4x10 cm measured h ere can be used together with 



the peak H2 column density measurement of iMeier et al.l (120021 ) to make a better estimate the 
extinction, Av{HI + H2), toward the dust lane, as described above. The column density of 
molecular hydrogen ca n be estimated fro m the integrated ^^CQ emission, /r^o, assuming that 



19871) 



a£m-i)/Icoii-o) 



1 (IMeier et al.ll200ll ). and adopting (IBloemen et al.lll986l : IScoville et al. 



NH^icm-^) ~ 3 X 10^%o(^ km s'^). 



(5) 



In terms of the CO flux, Fco, 



the direction of NGC 5253, the radial velocity in the heliocentric velocity frame is nearly identical 
(within 1 km s^^) to the L ocal Standard of R est frame, so we do not explicitly transform our 21-cm data 



into the LSR frame used by Meier et al. ( 20021 ). 
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NH,{cm-'^) c^Sx IO^OG {K Jy-^) Fco{Jy km 3-^] 



(6) 



where G is the effective gain of the a perture synt l iesis b eam, approximately 0.46 K Jy~^ 
for the beamsize and frequency of the iMeier et al.l (120021 ) Owens Valley Radio Observatory 
{OVRO) data at 230 G Hz. The peak CO f luxes for the four molecular clouds to the east of 
the nucleus measured by lMeier et al.l ( 120021 ) range from 6.9 - 16.9 Jy km s~^, corresponding 
to molecular hydrogen column densities = 0.95 - 2.33 x 10^^ cm~^. This may be a 
lower limit to the true molecular column density if the CO-to-//2 conversion factor in the 
low-r netallicity gas of NGC 5253 is larger than in the Milky Way, i.e., Xnq c 5253/-^Miy > 1 
(e.g., iMaloney fc Blacklll988l : IVerter fc Hodgdll995l : llsrael fc Burton! Il986l ). The maximum 
implied visual extinctions, based on peak ^8^^ beam-averaged H I and H2 column densities 
are Ay = 4.4 - 5.9. However, as noted by lCalzetti et al.l (119971 ) and evidenced by Figure [TT| 
the extinction derived from optical emission line ratio reddening maps is highly variable over 
small scales. These values range from Ay — 0.5 mag to Ay = 2 mag toward most of the star 
clusters. Ay may be > 9 mag toward the central embedded cluster, consi stent with estimate s 
of 8 mag < Ay < 25 mag from the 9.7 fim silicate feature optical depth (lAitken et al.lll982l ). 
Because the H I and CO interferometer measurements provide only beam-averaged extinction 
estimates over ~ 100 pc scales, localized regions may exhibit larger extinctions than the 4.4 
- 5.9 mag obtained from this approach. 



5.2. Comparison with the Hot Ionized Gas 

Figure [T2] is a composite 3-color plus contour image of NGC 5253 depicting the optical 
6450 A stellar continuum (blue), 0.3-8 keV X-rays (green). Ha (red) and t otal H I column 



density (contours, as in Figure [2]). The X-ray data, previously published in [Summers et al. 



(120041 ). was retrieved from the Chandra X-ray Observatory aichiYes (proposal number 02600546) 
and adaptively smoothed to a minimum signal-to- noise ratio of 2 for presentation. Figure [12] 
reveals the relative morphologies of the stars and the cold, warm ionized, and hot phases 
of the ISM in NGC 5253. The X-ray, Ha, and H I emission extend preferentially to the 
south and west of the nucleus, suggesting expansion of starburst heated material along these 
directions. If the H I plume is inflowing or interacting from the southeast, it may inhibit 
expansion in that direction. 
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6. Conclusions 



Radio- wave aperture syntliesis observations in the 21-cm line of neutral hydrogen at 8" 
resolution have revealed a complex non-axisymmetric gas distribution and velocity field in 
NGC 5253. The total neutral atomic content of the galaxy is 1.4 x 10® Mq from single- 
dish observations for a distance of 3.8 Mpc. The neutral medium in NGC 5253 may be 
characterized as consisting of three components: a component containing the majority of 
the neutral medium distributed along the stellar distribution and showing evidence for low- 
amplitude rotation at the level of 15-20 km s~^; a redshifted H I plume extending to the 
south-southeast along the minor axis containing 20-30% of the H I; and kinematically cold 
H I shells or filaments extending along the minor axes and to the southwest, containing 
~10% of the neutral gas mass. Taken together, these different components result in very 
complex velocity structures. Here we will try to summarize out best understanding of these 
structures, starting with the relatively enigmatic H I plume. 



6.1. The Nature of the H I Plume 

The kinematically distinct H I plume consists of a significant fraction of the total H I 
emission from NGC 5253. In this paper we have discussed three possible origins for the H I 
plume: outfiowing gas, infiowing gas, and a separate H I cloud which is interacting with 
NGC 5253. 

The possibility of outfiow was motivated by similarity with velocity fields observed in 
other dwarf starburst galaxies showing H I outflow (e.g., NGC 1569, NGC 625, NGC 1705). 
However, the large fraction of H I (?»25%) in the anamalous velocity gas in NGC 5253 
is exceptional. Furthermore, strong outflows are usually accompanied by Ha and X-ray 
emission which is not observed at this location. Such a massive outflow conflned to a small 
solid angle without corresponding signatures in the ionized gas seems unlikely. 

The possibility o f inflow was motiv ated by the radial velocities of CO clouds within the 



minor- axis dust lane (iMeier et al.ll2002l ). Our high resolution H I observations indicate that 
the plume is a kinematically distinct feature from the H I associated with the molecular gas. 
The most compelling argument for inflow is the near coincidence in velocity of the H I plume 
with the galaxy's systemic velocity in the region where the two are spatially coincident. If 
the H I plume is infalling, it must be from the far side of NGC 5253 given the plume's 
velocity fleld (i.e., the most redshifted gas lies farthest from the galaxy), consistent with the 
lack of extinction observed in this vicinity despite the high H I column density. 



Finally, we have the possibility that the H I plume is a distinct structure, perhaps a 
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remnant from a recent interaction, seen in projection against the main body of NGC 5253. 
In this scenario, a recent gravitational encounter, perhaps with another dwarf galaxy in the 
same group, produced the H I plume and provided the trigger for the current extreme burst 
of star formation. Other dwarf galaxies, both bursting and non-bursting are accompanied 
by pure H I structures. Given the relatively large gas mass associated with the H I plume, 
this would appear to be the most likely explanation. 



6.2. A Comprehensive View of the H I in NGC 5253 



Our high resolution view of NGC 5253 has shown us that NGC 5253 has many simi- 
larities to the best studied dwarf starburst galaxies. H I concentrations near the minor axis 
dust lane are coincident in loca tion and velocity with at least two of the molecular clouds 



seen in ^^CO (IMeier et al.ll2002l ). The H I extension to the south-southeast, designated here 



the H I plume, appears to be spatially and kinematically distinct from the gaseous medium 
surrounding the dust lane molecular clouds along the southeast minor axis. We present evi- 
dence that the redshifted H I plume is likely to be either a dynamically distinct tidal remnant 
or an inflow on the far side of the galaxy. H I extensions to the east, west, and southwest 
show remarkably similar morphologies as the Ha shells and filaments. In particular, there is 
a close m orpholo g ical s imilarity between the "quiescent" Ha shell to the west of the nucleus 
noted by iMartiru (119981 ) and an 800 pc diameter shell-like feature seen in H I. This may be 
an example of a starburst-powered shell stalled by the mass of surrounding neutral gas. 

These complex H I features and kinematic signatures, when considered together with 
distributions of the molecular, ionized, and hot coronal media, are consistent with a star- 
bursting low-mass galaxy in the initial phases of generating a galactic wind. Unlike more 
massi ve galaxies such as M 82 or NGC 3079 which sh ow monolithic, well-coUimated outflows 
(e.g.. [Strickland fc HeckmanI 120071 : ICecil et al.ll2002l . and references therein) several non-co- 
spatial star formation events over the las t ^ 20 Myr appear to be required to explain the 



multiple outflows observe d in NGC 5253 fjSummers et al.ll2004j : [Strickland fc Stevens! Il999 



Caldwell fc Phillipslll989h . The initiating and sustaining cause of the recent star formation 
activity may be an inflow or interaction of atomic gas along the "plume", identifled here 
for the flrst time. As such, NGC 5253 joins a long list of starburst galaxies where recent 
gravitational interactions are implicated as the trigger for extreme star formation bursts. 
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Fig. 1. — 21-cm emission line profiles (top) from the Parkes 64 m telescope (dotted line), the 
VLA D configuration data (solid line) and the VLA B+C configuration data cube (dashed 
line). The lower panel shows the 21-cm absorption line profile from the longest VLA baselines 
against the con apact thermal radio continuum source associated with the super star cluster 
in the nucleus ( iTurner et al.lll998l ). 
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Fig. 2. — NGC 5253 H I total intensity map (contours) with optical 6450 A continuum image 
(greyscale). Contours levels show H I beam-averaged column densities of 6, 12, 18, 24, 36, 
52, and 64 xlO^° cm~^. An 8" diameter (148 pc) cross marks the location of the dominant 
compact central star cluster and radio "supernebula" (ITurner et al.lll998l ). Note that the 
contours at the positio n of this feature ind icate absorption. Circles mark the positions of 
CO clouds reported by iMeier et al.l (120021 ). The inset at lower left shows the size of the 
9.0" X 7.6" synthesized VLA beam. 
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Fig. 3. — H I column density of NGC 5253 (contours) as in Figure [2] overlaid on an Ha 
image (greyscale). Note the correspondence between the Ha shell and H I arcs to the west 
of the nucleus. 
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Fig. 4. — H I column density of NGC 5253 (contours and greyscale) at a lower resolution than 
in Figure [2] in order to highlight the lower column density features. Contours correspond 
to beam-averaged H I column densities of 2, 4, 8, 16, and 32 xlO^° cm~^. Note the H I 
extensions to the east and west that coincide with Ha filaments. H I extensions along the 
major axis to the southwest also have Ha counterparts. 
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Fig. 5. — H I 21-cm maps of NGC 5253 in every 2.58 km velocity channel containing 
emission. The greyscale range shows beam-averaged column densities of 1.0x10^° cm~^ 
(white), equivalent to 2x the rms noise, to 5.0x10^'' cm~^ (black). The cross marks the 
location of the central super star cl uster and th e ellip ses are the molecular clouds A through 
E (from left to right) identified by iMeier et al.l (120021 ). 
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Fig. 5. — continued 
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Fig. 6. — H I column density channel maps of NGC 5253 (contours) overlaid on an Ha image 
(greyscale). The contours show fluxes of 2.4, 3.6, 6.0, and 8.4 xlO^^ Jy beam~^, equivalent 
to beam-averaged H I column densities of 1.0, 1.5, 2.5, and 3.5 xlO^° cm~^. 
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Fig. 7. — A 21-cm intensity-weighted velocity map of NGC 5253 (i.e., 1st moment; contours) 
overlaid on an optical 6450 A image (greyscale) . Heliocentric isovelocity contours are labeled 
in km s~^. Note the large velocity gradient in the southeast of the galaxy. Note also the lack 
of a strong gradient along the major axis of the galaxy. 
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Fig. 8. — A 21-cm intensity-weighted velocity map of NGC 5253 (i.e., 1st moment; contours) 
overlaid on a 21-cm total intensity image (0th moment; greyscale). Heliocentric isovelocity 
contours are labeled in km s"^. The scale bar shows H I surface brightness in units of Jy 
beam~^ m s~^. 
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Fig. 9. — A position- velocity diagram taken along the major axis at position angle 50°. Note 
the overall trend from left to right of more positive velocities, and the large velocity range 
near the center of the galaxy. 
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Fig. 10. — A position-velocity diagram taken along the minor axis at position angle 140°. 
Note the large, continuous, asymmetric velocity range near the center of the galaxy. 
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Fig. 11. — Extinction, Ay, as determined from the H I colum n density (cont o urs) a nd the 
Balmer emission-hne decrement (greyscale with colorbar; from ICalzetti et aLl (119971 )). The 
labeled contours correspond to Av{HI) =4, 3.5, 3, 2.5, 2, 1.5, 1 and 0.5 mag assuming a 
Galactic gas-to-dust ratio. The extended region of high extinction in the dust lane southeast 
of the nucleus has Av{Ha/ H (3) ~ 1.5 and is coincident with the second highest peak in the 
H I column density at Ay{HI) ~ 3.5. The H I column density is unreliable in the vicinity 
of the central radio supernebula (marked by a cross) owing to H I absorption. 
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Fig. 12. — A 3-color plus contour view of NGC 5253 showing the stars (green), Ha (red), 
0.3-8 keV X-rays (blue) and total H I column density (contours, as in Figure [2]). 
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